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1. Introduction 

The problem of the duration of the GRB prompt phase and its connection with their possible 
origin(s) is almost as old as their discovery. With the advent of Swift and the possibility to detect 
the prompt emission also in X-rays for some long GRBs, the burst durations extend over 6 orders 
of magnitude. This is only one of the multiple facet of then - diversity as their variability extend 
over 6 decades and the energy over 8 decades if we consider the afterglow phase. Based on their 
duration they have been categorized in 2 classes [1,2], respectively the short and long GRBs (here¬ 
after sGRBs and lGRBs). Though this classification is based on both their spectral and temporal 
properties, it uses only properties in the observer frame, being a problem for objects that have been 
observed from redshifts 0.1 to almost 9 (see also [3, 4] for an analysis of sGRBs in the rest frame). 
However, this categorization has proven to be useful, and lead to strong evidences for an association 
of sGRBs with compact binary coalescences [5], and lGRBs with the collapse of a very massive 
star to a black hole [ 6 ]. This classification is also correlated with studies on the localization of the 
events with respect to the host galaxy ( e.g. [7]). 

The collapsar model can explain the amount of energy needed to supply a 1GRB, and it is quite 
effective to explain many of the properties of these sources, such as the supernova associations 
[ 8 , 9] or the observation of stellar winds around the source [10, 1 1]. However, the lack of hydrogen 
in the envelope lead to the proposal that the progenitor of lGRBs should be a Wolf-Rayet star (e.e. 
[ 12 ]). 

This explanation is much more difficult with the so called ULGRB s (Ultra-Long Gamma-ray 
Bursts) as their duration is much larger than 1000s, and even above 20000s. To date we know 4 
candidate ULGRBs, namely GRB 110209A [13, 14, 15], GRB 101225A [15, 16], GRB 121027A 
[15, 17] and GRB 130925A [18]. 

The question of the origin of these events has been much debated during this conference and 
in the literature. Gendre et al. [13] have proposed that a blue supergiant as the progenitor could 
provide the matter reservoir needed to feed sources of such extreme duration. Levan et al. [15] 
reached the same conclusion, while Margutti et al. proposed a new class of soft and ultra-long 
events [19]. Recently [20] made a careful statistical analysis of set of long GRBs well observed 
over very long time and found that the probability that ULGRB s are an extension of the 1GRB 
class is very low, therefore they conclude that they form a different population. However, Virgili 
et al. [21] reached the opposite conclusion, and claimed that ULGRBs are the tail of the 1GRB 
distribution rather than event with different progenitors. 

One can also try to use other estimators of the GRB duration, such as modeling the X-ray light 
curve and using a theoretical model to estimate the duration of the central engine activity. The 
result of [22] is that the activity of GRBs span the range between 0.1 and 10 6 s, rather the observed 
peak distribution around 30s observed using the usual T 90 measure. 

Another point to add to the analysis is that all the four ULGRBs observed to date display a 
thermal emission during the whole prompt phase, a very unusual feature in GRB spectra. 

Whatever the case, ULGRBs being very long GRBs or a different classes, we are left with 
several questions. If ULGRBs are very long lGRBs, then the need arises to explain with the same 
mechanism bursts that last few seconds to more than 20000 s and with thermal emission in few 
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Figure 1: Template light curve with the notation in use in [20]. We present two examples. Above: one single 
break. Our method allow to discriminate if this break is prompt related or jet related. Below: a complex 
3-break light curve. In this example, the start of the steep decay is the second segment if a .2 > 2.2 (see more 
details in [20]). 


(usually very long duration) cases. If ULGRBs form a different class, then a credible new progen¬ 
itor, and/or physical mechanism has to be found to explain it. 

Here we present the results of a statistical analysis aimed at determined whether or not UL¬ 
GRBs, as a whole, form a different population or not. 

2. The new duration measure T x 

Following the work done in [20] we use the online Swift XRT GRB catalog 1 [23]. To model 
the overall shape of the light curves multiple power law segments /(f) = kt 01 are assumed. Flare 
episodes are considered as extra components in the XRT GRB catalog analysis, and removed during 
the estimation of the power law parameters. We have automated the analysis to extract the start time 

1 http://www. swift.ac ,uk/xrt_live_cat/ 
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Data 


Figure 2: The probability density function of the logarithm of Tx from [20]. We have superimposed three 
different fits; red curve: log-normal model on the whole sample; blue curve: log-normal distribution, sample 
restricted to Tx < 300s; brown curve: generalized extreme value distribution, whole sample. 


of the steep decay that we define to have a > 2.2. The advantages of this method are 1) that it uses 
only data from a single instrument, therefore it is spectrally consistent, 2 ) that it enables us to 
study GRB of duration larger than 1000s over several orbits. We have removed GRBs for which 
the follow-up starts more than 500s after trigger time, and we do not consider flares or late-time 
activity. The method is illustrated in Fig. 1 and in [20]. 

We are left with a final sample of 207 events, out of the 243 taken from the catalog. 

We have a specific case, GRB 130925A that has a claimed T 90 duration of more than 20ksec, 
but our Tx gives 149sec. This is probably due to a strong flaring activity. In any case we have 
considered for the analysis the two cases, with and without GRB 130925A. 

The mean duration of Tx = 137s and the median 119s. The minimum duration is 49s (but 
depends on the construction of the sample), while the maximum is 25400s. We have tested the 
different distributions: 

• Log normal distribution 

• Log normal distribution excluding data above 300s 

• Generalized Extreme Value Distribution 

The distribution of the logarithm of Tx as well as the above mentioned fits are plotted in Fig.2. 

The log-normal distribution on the whole set gives a very bad fit, and the null hypothesis is 
rejected to a high level (its probability is 2 x 10 5 ). If we truncate the distribution and retain only 
the events above 300s, we do not get an acceptable fit either. Interestingly, it is by removing the 
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events above 300s that we fit nicely the distribution. For example, keeping only events for which 
T x < 300 s, we get % 2 = 13 for 15 d.o.f. (while for T x < 2000s we got % 2 = 150 for 14 d.o.f.). This 
could indicate that our limit for ULGRBs (10 4 s) is too conservative, and a value of Tx > 10 3 s 
should be more representative of this new class. 

We get a better tit than the log-normal law with the GEV distribution. Flowever it is still non 
ideal because of an excess of events with Tx around 400-500s. The parameters are /./ = 1.99 ± 0.02 
(the location, i.e. 97.7s), the scale a = 0.19 ±0.01, and the form factor t, =0.17 ±0.05. Using 
this probability we get a prediction of 21 events above 400s, while 19 are observed. Flowever, the 
same law gives the probability to get a point above 10 4 s to be 2 x 10 3 . 

We tested the addition of GRB 130925A by adding its claimed duration of 20ks [18] obtaining 
a very similar result. 

In conclusion these results indicate clearly that the presence of one (or two) GRB above lOksec 
is an outlier. 

3. Discussion 

In the previous section we show that the light curve of ULGRBs is markedly different from 
that of lGRBs. Not only they display burst activity for very long times, but also there is a "regular" 
afterglow decay after times larger than lOksec (see. e.g. the statistical analysis performed in [18]). 

Another set of evidence come from the case of GRB 130925A that showed constant radio 
emission over four months, the clear sign of a wind environment. Finally, as already noted, ah 
ULGRBs display a thermal component. As an example, GRB 130925A has been observed up to 
10 7 s by XMM-Newton, with a IkeV thermal spectrum, a temperature slowly decreasing, which led 
[18] to conclude that it could be the sign of a hot cocoon. 

The case of GRB 130925A rises the question of whether or not we should have considered 
late flares in the burst duration. Late flares can be the signature of a sustained activity of the central 
engine. On the other hand, they can be the signs of the injection of slow blobs in the fireball that 
appear at later times. In any case, whatever the conclusion, the activity of the central engine for 
more than 20ksec needs an explanation. 

Another question that arises is that of the progenitor: should we invoke a different source of 
energy? Our work show that ULGRBs are not an extension of normal 1GRBS, and that they display 
also a thermail component. Accretion from an extended source seems a natural explanation, as 
already noted [13, 24]. 

4. Conclusions 

We have presented evidences that ULGRBs form a class different of that of "regular" lGRBs. 
Accretion from an extended source is a natural explanation, as already noted, though probably not 
the only one. It is tempting to see ULGRBs as large, extended very low metallicity stars, fossils of 
PopIII stars, e.g. blue supergiants. 

The fact that the detected ULGRBs are ah relatively nearby is not really suprising, as they 
can’t be detected at redshift larger than 1.5 in the very best case as shown in Lig.3. If we want 
to detect more of these guys, new experiments, such as the SVOM satellite [25], should use new 
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Figure 3: Detection significance expressed in signal to noise ratio units for a GRB 111209A like event as a 
function of its redshift. We present the detection regime of such bursts using light curve bin sizes of 1 minute 
(purple line) and 2 minutes (blue line); and a burst 10 times fainter with 2 minutes bin size (black line). The 
dotted red line indicates our detection threshold (5 sigma level), showing that even in the optimal case the 
maximum redshift for detection is z = 1.4-1.5. From [13]. 


trigger criteria that maximize the probability to detect distant events. As ULGRB s might well be 
the swansong of PopIII supergiant stars, this can be considered as a priority. 
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